Double butterfly spectrum for two interacting particles in the Harper model 
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We study the effect of interparticle interaction U on the spectrum of the Harper model and show 
that it leads to a pure-point component arising from the multifractal spectrum of non interacting 
problem. Our numerical studies allow to understand the global structure of the spectrum. Analyt- 
ical approach developed permits to understand the origin of localized states in the limit of strong 
interaction U and fine spectral structure for small U. 
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Recently a great deal of attention has been devoted to 
the investigation of incommensurate systems exhibiting 
singular continuous spectrum with many interesting mul- 
tifractal properties (see e.g. Among the physical 
models, one of the most popular is the Harper model of 
electrons on a two-dimensional square lattice in the pres- 
ence of a perpendicular magnetic field [|,||] . This system 
can be reduced to the study of a rather simple model 
of particle dynamics on a one-dimensional quasiperiodic 
lattice. The energy spectrum exhibits multifractal prop- 
erties and the band spectrum for rational values of mag- 
netic flux looks like a butterfly. In spite of the academic 
character of such a model, experiments have been per- 
formed during the last ten years exhibiting this multi- 
fractal butterfly structure. One of the first among them 
has been performed in 1985 using superconducting net- 
works H and more recently experiments with superlat- 
tices also allowed to observe the fist hierarchical steps of 
multifractal butterfly structure j?J . 

The deep understanding of such an intricate spec- 
tral structure attracted interest of mathematicians and 
physicists who developed new approaches for its inves- 
tigation such as non commutative geometry 0, pseudo- 
differential operators M , functional analysis |1C| , renor- 
malization group approach 1 1 1 13] , thermodynamical for- 



malism 1 13 . All these tools allowed to study the problem 
on rigorous mathematical ground and to understand the 
properties of eigenstates. For example using the dual- 
ity between momentum and spatial coordinate Jl4[| , it is 
possible to prove rigorously the existence of localized or 
delocalized states |T^,|l6|. It was also found that quan- 
tum systems which are chaotic in the classical limit may 
have quite unusual properties in the presence of underly- 
ing quasiperiodic structure p7|,p]p^] . 

All the works mentioned above were done for one par- 
ticle dynamics. However even from the physics of the 
original Harper model, it is clear that the interaction be- 
tween electrons on the square lattice in the presence of 
magnetic flux plays an important role. Therefore it is 
natural to address the question of the influence of inter- 



action on multifractal spectrum. The most simple exam- 
ple of such a case is an interaction between two particles. 
Recently it has been found that in the case of random 
potential even such simple model has a number of unex- 
pected properties Ga|. For example repulsive/attractive 
short range interaction leads to appearance of effective 
pair states in which two particles propagate together on 
a distance much larger than the one-particle localization 
length without interaction. Surprinsingly the first nu- 
merical studies of interaction effect in a quasiperiodic 
potential showed an opposite tendency fUjfl . Namely, re- 
pulsive/attractive interaction leads to the appearance of 
localized states while in the absence of interaction mul- 
tifractal spectrum generated quasidiffusive spreading of 
wave packets on the lattice. However, the numerical 
approach used in jl9| allowed to study only the wave 
packet evolution while the structure of the spectrum it- 
self was not directly accessible. Therefore to understand 
the spectral structure and the nature of eigenstates we 
performed numerical simulations by direct diagonalisa- 
tion based upon Lanczos algorithm. 

As a basic model for our investigations we consider the 
model of two interacting particles (TIP) in the Harper 
problem described by the following eigenvalues equation 

(2Acos(7ni + 2Acos(7n 2 + (3 2 ) + U5 nun2 )<p nun2 + 

<Pni + l,n 2 + ( y 5 ni-l,ri2 + ( Pn 1 ,n 2 + l + 'Pni^-l = -EVni,n 2 

(1) 

where the parameter 7 characterizes the quasiperiodic 
lattice for the one-particle problem. Without interaction, 
each particle moves in quasiperiodic Harper potential and 
7/27T = (f)/(t>o = a is the ratio between the magnetic flux 
within one unit cell of the square lattice and the flux 
quantum </) = hje. The parameter a plays the role of 
an effective Planck's constant so that a i— > corresponds 
to the semiclassical limit. The two parameters /3i.2 are 
related to the quasimomentum components in the non 
interacting problem. The parameter A characterizes the 
strength of the quasiperiodic potential and for the case 
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of electrons on a square lattice A=l |. However from 
mathematical point of view it is also interesting to study 
the different regimes with A < 1 and A > 1. Strong 
analytical and numerical evidence has been given that 
the spectrum is pure point and the states are localized 
when A > 1 while for A < 1 the spectrum is continuous 
with extended eigenstates 14 10|Jl],po|. The strength of 
the short range on-site interaction is characterized by U. 
We concentrate our investigations on the case A = 1, 
/?i.2 = P when for U = the spectrum is multifractal for 
irrational values of 7/277. We consider only the part of 
the spectrum corresponding to the symmetric TIP states 
since antisymmetric configuration is not affected by on- 
site interaction. 

In the absence of interaction, the corresponding two 
particle spectrum results of the superposition of two one- 
particle spectra of the Harper model and is shown in Fig. 
1 (a). Comparing with the one-particle spectrum (Hof- 
stadter's butterfly), we can remark that the spectrum 
becomes much more dense near the centers of the bands 
and subbands but still the gaps in the spectrum survive 
on all energy scales. 



When increasing the strength of the interaction U, the 
spectrum is splitted into two butterflies which are slightly 
shifted one respect to the other. However one of them 
remains almost at the same place corresponding to the 
non interacting case of Fig. 1 (a). The shifted butterfly 
moves to the right since the repulsive interaction U > 
gives global increase of energy. A typical case U = 1 of 
double butterfly spectrum is presented in Fig. 1 (b). 

The main features which can be immediately observed 
in this figure are the smoothness of the edge of the shifted 
butterfly, the less dense character of its spectrum and 
the filling of some internal energy gaps (sec for example 
near a — 0.6 and E = —1.5). However, the gaps in the 
spectrum still exist on all scales. 

The shift of one butterfly and almost unchanged form 
for the other at moderate values of interaction U can 
be understood in the following simple way. For that we 
choose small values of flux a < 1 and use the pertur- 
bation theory in U on the basis of harmonic oscillator 
functions to get analytical expressions for the Landau 
sublevels at the spectrum edge. Without interaction, the 
band edge is given by E±(a) = ±8=F47ra(m 1 +m 2 + 1) ± 
ir 2 a 2 (2+ (2mi + l) 2 + (2m 2 + l) 2 ) /4 + 0(a 3 ), which is 
superposition of two Hofstadter butterflies in semiclas- 
sical regime pj] . The integers m\,m,2 are the Landau 
quantum numbers for oscillator states near the bottoms 
of potential minima. If two particles are located in dif- 
ferent minima, the interaction between them is negligibly 
small and the energy levels are not shifted by U. These 
energy states correspond to non shifted butterfly with 
dense spectrum since there are many states when TIP 
are separated from each other. If TIP are located in 
the same potential minimum, the interaction gives en- 
ergy shift which in the first order of perturbation theory 
is AE± = U^fa for 7711,2 = and 7711,2 = (0; 1) being in 
good agreement with numerical data for U < 1 as can 
be seen on Fig. 2 (a). This shows that the shifted but- 
terfly corresponds to the case when the two particles are 
located near each other. The density of such states is 
smaller than in the case when particles are far from each 
other and that is why the shifted butterfly is less dense. 





Fig. 1 : Spectrum of two particle Harper problem (a : 
up), with U = obtained for rational values 0/7/277 = 
a = p/q with q < 19; (b : down) with U — 1 and q < 23. 



Fig. 2 : Energy band edges (a) U—0.4, dots are numer- 
ical data and solid curves are perturbation theory results 
(see text); (b) U=10, dots are data from Fig. 4 and solid 
curves are given by theory described in the text. 

Direct analysis of eigenstates for irrational flux values 
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(which are approximated by a continuous fraction ex- 
pansion) also shows that the states in the shifted part 
correspond to the situation where two particles stay near 
each other. However, contrary to the TIP in a random 
potential, the particles here cannot propagate together 
and stay exponentially localized near the origin as it can 
be seen with the typical 3-D plot of Fig. 3. 



two times smaller than the original spectrum at U = 0. 
The center of this band is located at energy E « U. The 
typical example of global spectrum is shown in Fig. 5. 



Fig. 3 : Semilog plot ofW ni ,n 2 — l^ni.nal f or localized 
(E = -1.3376,-10 < ]nW < = 5.9,£ = 193 (a 

: left)) and delocalized (E = -1.7368,-10 < h\W < 
— 3,£ = 214, £o = 12.5 (b : right)) eigenstates at U = 
l,a= 34/55, 13 = y/2. 

We also investigated the structure of eigenstates in 
the more dense part of the spectrum (non shifted but- 
terfly). In that case the eigenstates are delocalized and 
quite similar to those corresponding to the non interact- 
ing case. Here the two particles mainly spread quasidif- 
fusively along the quasiperiodic lattice and interaction is 
not important for them. This structure of localized and 
delocalized eigenstates is in agreement with the numeri- 
cal study of wave packet dynamics performed in Jl9| . 

The properties of eigenstates can also be analyzed 
with the help of the inverse participation ratio (IPR) 
£ = (E n ^rei,tu) _1 - value for different energies 
is shown in Fig. 4 for a — 34/55. In agreement with the 
above discussion, the localized states with small £ corre- 
spond to the part of the shifted butterfly with less dense 
spectrum while the unshifted butterfly is associated to 
large £ with delocalized states. It is interesting to deter- 
mine the IPR £o in the non interacting eigenstates basis. 
Such approach has been quite useful for TIP in a random 
potential |^2| . It is interesting to note that the situation 
for TIP in the Harper model is quite different. Namely, 
the delocalized states have very small value of £o while 
the localized ones are delocalized in the non interacting 
eigenstates basis and have very large £o (see Fig. 3). This 
result once more shows that delocalized states correspond 
to almost non interacting particle propagation while lo- 
calized states appear only due to interaction which can 
be even repulsive (Fig. 3a). 

With further increase of U the shifted butterfly goes 
on moving to the right and becomes more and more de- 
formed. Starting from interaction strength U > 10, this 
butterfly is transformed into a spectral band with width 



Fig. 4 : Inverse participation ratios £ vs eigenenergies 
E shown at £ = 2; U = 1, a = 34/55, < (i < 2tt. 



Fig. 5 : Same as in Fig. 1 with U = 10 and q < 28. 

The physical reason for the appearance of such sepa- 
rated spectral band can be understood in the following 
way. For strong U, there are states for which TIP are 
localized on the same site so that = n. According to 
(1) the energy of the states is E n = 4Acos(7n + (3) + U. 
The transition between these states can be obtained with 
first order perturbation theory in 1/U which gives the ef- 
fective eigenvalue equation : 



(4A cos(7n + {3) + U) (f>„ + V e 



eff ! 



Jn+l 



+ (f>n-l) = E(f> r , 



(2) 



Here V c s is the hopping between such states due to vir- 
tual transitions via states with n\ ^ U2- For U 3> 1, 
the energy difference between diagonal and off-diagonal 
states is very large and therefore V e g ~ 1/U. The equa- 
tion for diagonal eigenstates has the form of Harper equa- 
tion with A e ff = 2A/V e g ^> 1. Due to that these states 
are exponentially localized so that particles stay near the 
origin. In some sense, the interaction renormalizes the 
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constant A i— > A G ff in the Harper equation for pair of par- 
ticles. For strong U, the renormalized A c ff is much larger 
than 1 that, according to the Aubry duality |14|, leads to 
localization of TIP pairs in quasiperiodic potential. Our 
conjecture is that in a sense A e g remains larger than 1 
even for moderate values of U ~ 1. In a sense interaction 
breaks Aubry duality leading to appearence of localized 
TIP phase. However more rigorous analytical confirma- 
tions of this conjecture are desirable especially keeping 
in mind that in a random potential the interaction with 
U ~ 1 leads to derealization of TIP pair states. The 
accurate expressions for the TIP energy edges of shifted 
spectral band can be found using semiclassical analysis at 
small flux values by methods developed in [p| j23[ . The de- 
tails of computations will be given elsewhere ]p4[ . For the 
case of Fig. 2 (b), they give E = 6.0 + 0.59*27ra + O(a 2 ) 
that is in good agreement with numerical data (Fig. 2 
(b))- 

For the part of the spectrum represented by un- 
shifted butterfly at U ^S> 1, the eigenstates become more 
and more similar to asymmetric TIP configuration i.e. 

= signK - n 2 ) ( X M - X £W) /V2, where 

X~s are one-particle eigenf unctions. Due to that, the ef- 
fective interaction becomes quite small and the unshiftcd 
butterfly at large U (TIP are in different wells) looks 
very similar to the one at U = 0. The main difference 
is the splitting of Landau sublevels which appears due 
to effective small interaction between particles located 
in the same well. According to the expression for 4> ni ,n 2 i 
such splitting can take place only when Landau quantum 
numbers are different (mi 7^ 7712) so that % 7^ X • As 
the result the first sublevel with mi. 2 = is not splitted. 
For non interacting part, the edges are given by the same 
E±{a) as for U — (see above) while for interacting case, 
the additional shift is 6E(a) = -8?ra/(Z7 + 4) (see [||). 
These analytical expressions are in good agreement with 
numerical results as shown in Fig. 2 (b). 

In summary, 20 years after 0] our investigations of 
spectra and eigenstates for TIP in the Harper model (1) 
show that repulsive/attractive interaction leads to ap- 
pearance of localized states. Our conjecture is that due 
to Aubry duality breaking a localized TIP pair phase 
appears at arbitrary small interaction strength. At the 
same time we expect that this breaking is absent for TIP 
on the 2d-lattice with magnetic flux. However, the later 
model requires separate investigations [ p4f . 
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